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Abstract: The LHC accelerator complex will be upgraded between 2020-2022, to the High- 
Luminosity-LHC, to considerably increase statistics for the various physics analyses. To operate 
under these challenging new conditions, and maintain excellent performance in track reconstruction 
and vertex location, the ATLAS pixel detector must be substantially upgraded and a full replace¬ 
ment is expected. 

Processing techniques for novel pixel designs are optimised through characterisation of test struc¬ 
tures in a clean room and also through simulations with Technology Computer Aided Design 
(TCAD). A method to study non-perpendicular tracks through a pixel device is discussed. Com¬ 
parison of TCAD simulations with Secondary Ion Mass Spectrometry (SIMS) measurements to 
investigate the doping profile of structures and validate the simulation process is also presented. 
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1. Introduction 

The Large Hadron Collider (LHC) [0] will be upgraded to the High-Luminosity LHC (HL-LHC) 
during the phase II upgrade expected in ~2022. These modifications will result in an increase 
in occupancy and of radiation damage to the ATLAS [0] Inner Detector (ID), shown in figure 0. 
From the interaction point outwards, the current ATLAS ID consists of the pixel sub detector, the 
Semi-Conductor Tracker (SCT) and the Transition Radiation Tracker (TRT). In 2014 an extra layer 
called the Insertable B-Layer was installed between a smaller beam pipe and the first layer of the 
pixel sub-detector. During the phase II upgrade, the ATLAS experiment will replace the current 
ID, with an all silicon tracker composing of pixels and strips. This will result in a greater area 
covered by pixel detectors compared to the present design. The expected fluence for the inner-most 
pixel layer (at approximately 4mm) will be 2 x 10 16 n eq cm~ 2 , where 1 n eq is 1 MeV neutron 
equivalent fluence. Consequently the design requires an entirely new ATLAS pixel sub-detector 
with improved pixel devices which are radiation hard, have slimmer edges and a better granularity. 

Through new pixel designs, such as active edges and alternative bias rail geometries, minimi¬ 
sation of inactive regions and reductions in efficiency loss can be obtained, allowing sensors to 
be placed adjacent to each other. This layout is preferable to shingling as it reduces the material 
budget of the ID as well as cooling requirements and power consumption. 

This article presents work performed in the context of this upgrade programme including the 
study of methods to analyse non-perpendicular particle incidences through pixel prototype devices 
in a clean room environment, and also the analysis of doping profiles for specifically manufactured 
samples with parameters used for ATLAS pixel devices. 
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Figure 1: The ATLAS Inner Detector of Run I (Image courtesy of CERN). 

2. Laboratory Characterisation 

Previous work on characterisation studies at LAL, Orsay have resulted in the redesign of the guard 
ring structure of planar pixel sensors, reducing the inactive region at the edge of the device from 
1100 jUm to 200 jtmi [0]. 

Devices are tested in a clean room environment with the USBPix system ||]], which is com¬ 
posed of a muti-IO board and an adapter card designed specifically for the readout-cards known 
as FE-I3 or FE-I4. The setup requires scintillators with photomultiplier tubes, one above and one 
below the device under test, for triggering on particles traversing the device either from a radioac¬ 
tive source, or from cosmic muons. The analogue signal from the photomultiplier tube is then fed 
into the wavecatcher [ 0 ] device which outputs a digital signal trigger pulse from the coincidence 
of the two signals. The results shown below are a proof-of-principle test on a well known device 
to explore the possibility of using the setup with horizontally offset scintillators, as illustrated in 
figures ^a] and [ZF]. to study non-perpendicular cosmic muon tracks through the sample. 

The scintillators are approximately 26 mm by 26 mm. The first measurement taken was with 
a full overlap of both the scintillators with the sensor. This “full-overlap” case assumes that a high 
proportion of the tracks are perpendicular to the sensor surface, however it should be noted that an¬ 
gles of up to 40° from perpendicular are possible due to the size of the scintillators. Subsequently, 
the two scintillators were offset in the horizontal plane, in opposite directions, so that each scin¬ 
tillator overlapped with half of the sensor, but not with each other (this is the setup illustrated in 
the sketch shown in figure [2a| ). In this case, the angle of the track incidence on the sensor would 
be between 0° and 60°. The perpendicular tracks would be rare as and only pass through the cen¬ 
tre of the setup. Finally, the offset was increased to introduce a gap of 20 mm between the two 
scintillators (exactly the width of the sensor) to select tracks between approximately 34° and 61°. 

The mean cluster size, as shown in figure ^j, increases as the scintillator gap is widened hori¬ 
zontally, as expected. The advantage of this method is that it can be used to study non-perpendicular 
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Figure 2: (a) Sketch and (b) Photograph of the setup of offset scintillators, shown above and below 
the device being testes, to provide non-perpendicular tracks. 



Cluster 

Figure 3: Cluster distribution map for FE-I4 planar pixel sensor device for perpendicularly arranged 
scintillators (blue), half offset (green) and full offset (red). 


tracks through the sensor when test beams are not available. When combined with cluster size in¬ 
formation, this method could allow studies of the charge collection at various depths, and therefore 
various depletion regions, within the sensor bulk [0] . However, since there are fewer muons that 
pass non-perpendicularly and the trigger window is smaller, the method is time consuming, taking 
on the order of a full 24 hours to collect 1,000 events. Reconstruction of track paths is also not 
possible without a telescope. 
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3. SIMS and TCAD Simulation 


Motivation Results that have been obtained from test structure measurements can be used to 
develop reliable simulations of device, which, in turn, drive the development of new sensor layouts. 
This process is quicker and more cost effective than building multiple physical prototypes. It is 
vital, therefore, that the simulation is validated to be able to trust predictions. Studies of irradiated 
devices are highly important for the preparations of the ATLAS phase II upgrade and also in the 
context of the CERN RD50 programme [0], but the simulation must be first validated for non- 
irradiated samples. 

Secondary Ion Mass Spectrometry (SIMS) [0] is a process by which the impurities in the sur¬ 
face, and near surface region, of a sample can be measured [@|. These impurities are intentionally 
implanted into a silicon wafer to produce a doped region. For n-doping, impurities such as phos¬ 
phorus are used, where as boron is commonly used for p-type doping. The SIMS measurement is 
performed by sputtering a primary energetic ion beam onto the sample and measuring the produced 
ionised secondary particles via mass spectrometry. This process is destructive, leaving a crater in 
the sample, and therefore specific samples produced for this analysis are usually required. The ma¬ 
chine used for this study was the Cameca IMS 7F system, at the GEM AC laboratory in Versailles; 
further details can be found in reference fITTlIl . 

The SIMS process is shown in figure 0. The primary ion source comes either from a duo- 
plasmatron (1) or from a surface ionisation source (2). A primary beam mass filter is within the 
primary ion column (3), which leads to the secondary ion extraction transfer (4). The ion energy 
is analysed at (5), the mass at (6) and the secondary ion extractors are at (7) and (8). The count 
rate of the elements chosen for analysis is measured as a function of time, which is subsequently 
converted to depth profiles with the use of a mechanical stylus. For further details of the analysis 
process see reference [Ql . 

It should be noted that this method measures the total atomic dopant profiles only, which is 
the physical concentration of impurities in the sample; it does not give information about the active 
dopant profile, that which actively contributes to the electrical field within the sample. Therefore, 
the total dopant profile measured by the SIMS method is not expected to change after a device has 
been irradiated, but it is important that this is verified. 


3.1 SIMS Sample Preparation 

Samples are prepared by dicing uniformly n- or p-doped, non-thinned wafers of 6- or 8-inches with 
a diamond saw. The circular wafer is diced into four quadrants with four strips taken from the 
centre. Three quadrants are kept for the future in case further analysis should be performed. There 
is a layer of oxide present on the top of the sample; this layer is an insulator. The annealing process 
can sometimes cause an increase in the thickness of this oxide layer. Measurements are taken either 
through the oxide, or chemical oxide etching is performed to remove this layer. Samples for the 
study were manufactured at two facilities; n-in-n samples were produced at CiS, Germany [Ql 
and n-in-p samples were manufactured at VTT, Finland [Ql , both with a design oxide thickness of 
100 nm and 200 nm. High resistivity wafers corresponding to > 4 kil cm 1 , and low resistivity of 
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Figure 4: Diagram of the SIMS process. The details of each stage, indicated by the numbers, are 
described in the text. 

Table 1: Samples manufactured at CiS for analysis with the Secondary Ion Mass Spectrometry 
process. Various implantation doses and energies for the samples were chosen, with the aim of 
studying the total doping profile in the near surface region. 


n-in-n, CiS production, <100> orientation, thickness 380 gm 


Oxide thickness 


100 nm 


P implantation doses (atoms) 

10 13 cm" 2 

10 14 cm 2 

10 15 cm 2 

10 16 cm- 2 

Implantation energy (keV) 

130 240 

130 240 

130 240 

130 

Annealing 


4 hours, 975 °C 



0.25 £2 cm -1 , were selected. Four implantation doses ranging from 10 13 at. cm -2 to 10 16 at. cm~ 2 
were chosen for study, with two implantation energies, 130 keV and 240 keV. These values were 
chosen because they are close to the ATLAS pixel detectors fabrication parameters. Tables [T] and 0 
show the two production runs for the respective manufacturers for the 100 nm and low resistivity 
batches only. It should be noted that VTT used a different annealing process to CiS, however the 
details cannot be revealed due to confidentiality agreements. 

The samples with 200 nm oxide layer and high resistivity which were also produced for this 
analysis process, are for a further study and therefore will not be presented here. 

3.2 TCAD Simulation 

Simulation is performed with Technology Computer Aided Design (TCAD). Two versions have 
been used: Silvaco UHl and Synopsys [IT511. The Silvaco diffusion models considered are the 
following: Fermi model; 4 CPL model (clustering consideration); and PLS Solid State model. For 
Synopsys: Charged Pair Model; Charged Fermi model; Constant model; Charge React model. The 
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Table 2: Samples manufactured at VTT for analysis with the Secondary Ion Mass Spectrometry 
process. Various implantation doses and energies for the samples were chosen, with the aim of 
studying the total doping profile in the near surface region. 


n-in-p, AD VAC AM production, <100> orientation, thickness < 675 jum 


Oxide thickness 


100 nm 


P implantation doses (atoms) 

10 13 cm 2 

10 14 cm- 2 

10 15 cm 2 

10 16 cm 2 

Implantation energy (keV) 

130 240 

130 240 

130 240 

130 240 

Annealing 


3 hours, 1000 °C 



details for each model can be found in the respective user guides. The models use varying numbers 
and complexities of equations to compute the diffusion process within the simulation. Models 
are chosen for the appropriate level of accuracy and computational cost. As such, for Synopsys, 
the Constant model assumes a fixed diffusivity with no electric field and is therefore the least 
computationally expensive, while the Charge React model solves seven equations and has a higher 
computational cost. Simulations with Silvaco have been performed, but will not be presented here. 

The parameters for oxidation are included in the simulation IfTHl . instead of adding a layer of 
oxide on top manually, and the simulation steps are taken from the manufacturer’s own process. 

3.3 Results 

Simulation using the Synopsys Charged Pair diffusion model, and the SIMS data for the CiS low 
resistivity wafers batch have been compared, as can be seen in figure |5a|. There is generally a 
good agreement between them, however, a large discrepancy at the highest dose (10 16 at. cm—2) at 
greater depths has been observed. 

The same comparison was performed for VTT low resistivity wafers, shown in figure [5F]. Here 
the discrepancy for the highest dose is not observed. To study this, various diffusion models were 
used, such as the Charge React and the Constant model for the CiS sample in figures ^a| and |6Fi| 
respectively. The discrepancy for the two CiS doping profile samples of 10 15 at. cm~ 2 is understood 
to be due to a systematic issue, rather than an error with the simulation, but further study is required. 

Figures [7a] and [715| show the same comparison for VTT samples as described above for CiS, 
however results for the Charged React model are not shown due to similarities with the Charged Pair 
model and consequently the Charged Fermi model was studied. The comparison for the Charged 
Fermi model was poorer for the highest dose, and for the Constant model, all doses exhibited a 
significantly worse comparison. It was shown that for low and intermediate doses exhibit relatively 
good agreement, however for each model there is still a discrepancy for the dose of 10 16 at. cm -2 . 

As can clearly be seen, no model accurately simulates the SIMS data. Further work is required, 
possibly to merge two models together, combining the low depth and high depth simulations. 
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Figure 5: TCAD simulations performed for various implantation doses (at. cm~ 2 ) and energies 
(keV) for the Synopsys Charged Pair model, compared to SIMS measurements of (a) CiS and (b) 
VTT low resistivity samples. 




Depth [nm] Depth [nm] 

(a) (b) 

Figure 6: TCAD simulations performed for various implantation doses (at. cm~ 2 ) and energies 
(keV) for (a) Charged React and (b) Constant models, compared to SIMS measurements of CiS 
low resistivity samples. 


4. Summary and Conclusions 


Characterisation at LAL of novel pixel prototypes for the HL-LF1C, including a method to study 
non-perpendicular tracks through an FE-I4 sensor, has been discussed. 

Results have been presented for doping profile measurements performed with the SIMS method 
for samples from both CiS and VTT (ADVACAM). These measurements were compared with 
TCAD simulations for Synopsys and Silvaco. A discrepancy has been observed between simula¬ 
tion and data for the highest doping profile for all simulation models for CiS structures. Further 
studies for 200nm oxide layer and high resistivity wafers will be performed, and samples will also 
be irradiated to investigate what effect this will have on the total doping profile. 
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Depth [nm] Depth [nm] 

(a) (b) 

Figure 7: TCAD simulations performed for various implantation doses (at. cm~ 2 ) and energies 

(keV) for (a) Charged Fermi and (b) Constant models, compared to SIMS measurements of VTT 

low resistivity samples. 
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